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Abstract 

I evaluate the thermodynamic curvature R for fourteen pure fluids 
along their liquid- vapor coexistence curves, from the critical point to 
the triple point, using thermodynamic input from the NIST Chem- 
istry WebBook. In this broad overview, R is evaluated in both the 
coexisting liquid and vapor phases. R is an invariant whose magni- 
tude is a measure of the size of mesoscopic organized structiures in 
a fluid, and whose sign specifies whether intermolecular interactions 
are effectively attractive (R < 0) or repulsive (R > 0). I discuss five 
principles for R in pure fluids: 1) near the critical point, the attractive 
part of the interactions forms loose structures of size \R\ proportional 
to the correlation volume and sign of R negative, 2) in the vapor 
phase, there are instances of compact clusters of size \R\ formed by 
the attractive part of the interactions and prevented from collapse by 
the repulsive part of the interactions, and sign of R positive, 3) in the 
asymptotic critical point regime, the i?'s in the coexisting liquid and 
vapor phases are equal to each other, i.e., commensurate, 4) outside 
the asymptotic critical point regime incommensurate i?'s may be as- 
sociated with metastability, and 5) the compact liquid phase has \R\ 
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on the order of the volume of a molecule, with sign of R negative for a 
liquidlike state held together by attractive interactions and sign of R 
positive for a solidlike state held up by repulsive interactions. These 
considerations amplify and extend the application of thermodynamic 
curvature in pure fluids. 

Suggested PACS Numbers: 05.70.-a, 05.40.-a, 64.60.Bd 

1 INTRODUCTION 

Statistical mechanics and thermodynamics are based naturally in different 
domains. Statistical mechanics starts at the microscopic level with known 
intermolecular interaction potentials. Averaging with the Gibbs-Boltzmann 
distribution leads from this microscopic domain to the macroscopic domain 
where the laws of thermodynamics apply [H |2] . But there is an intermediate 
mesoscopic domain where essential fluctuation phenomena connected with 
phase transitions take place. The theme of this paper is getting information 
in this difficult domain using the thermodynamic curvature R. 

At least conceptually, statistical mechanics allows us to "build up" from 
the microscopic level to calculate properties in this mesoscopic regime. But 
this operation is generally difficult in practice. Less conceptually natural is 
thermodynamic fluctuation theory [H 12], in which we "build down" to the 
mesoscopic domain from the thermodynamic one. But information gets lost 
in the averaging yielding thermodynamics from statistical mechanics, and 
the idea of getting some of this information back seems counterintuitive at 
first. Here, I do this backtracking using R, which connects to intermolecular 
interactions O S] . 

In this paper, I calculate R along the liquid-vapor coexistence curve from 
the critical point to the triple point for fourteen pure fluids, in both the liquid 
and vapor phases. A schematic fluid phase diagram is shown in Figure 1, 
where (p, T) denote pressure and temperature, respectively, with subscripts 
"c" and "t" for critical and triple point properties, respectively. The ther- 
modynamic input for this broad survey calculation comes from the NIST 
Chemistry WebBook [5J, based on fits of experimental data. Although the 
contents of this data base will not always be optimal in any particular re- 
gion, or necessarily contain the latest experimental data, the NIST Chemistry 
WebBook represents the state of the art in representing data for many fluids 
over a large span of thermodynamic states. The calculations in this paper 
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Figure 1: A typical phase diagram for a material consisting of only one type 
of molecule. Three phases, solid, liquid, and vapor are separated by first- 
order phase transition lines. Of interest in this paper is the liquid-vapor 
coexistence curve connecting the triple point and the critical point. The 
asymptotic critical region is indicated with an open circle. 

mark the logical first step in understanding the broad behavior of R in pure 
fluids. 

This paper is arranged as follows. First, I summarize the method of 
calculation of i? in a number of thermodynamic coordinates. Second, I give 
an overview of the physical interpretation of R in pure fluids. This overview 
includes both what was known previously, and what was learned here. Third, 

1 give results for R calculated in fourteen pure fluids with the NIST Chemistry 
WebBook. Fourth, I have an Appendix with proofs of critical point properties 
of R. 

2 CALCULATION OF R 

In this section, I summarize how R is calculated in fluids in various coordinate 
systems. For a given thermodynamic state, R is an invariant, the same cal- 
culated in any thermodynamic coordinate system. However, an appropriate 
coordinate system can much simplify a calculation. 
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Consider a thermodynamic system consisting of one type of molecule, and 
with fundamental equation U = U{S, N, V), where U is the internal energy, 
S is the entropy, N is the number of molecules, and V is the volume [6]. 
Define as well the temperature, chemical potential, and pressure: {T, ^ 
{U^SyU^N, —Uy}, where the comma notation indicates differentiation. De- 
fined in Table 1 are the Helmholtz free energy A = A{T,N,V), the grand 
canonical potential fl(T, fi, V), and the Gibbs free energy G = G(T,p, N). 

The thermodynamic entropy information metric (A£)^ is defined in terms 
of the fluctuation probability of an open subsystem with fixed volume V of 
an infinite reservoir in a reference state [H |2]: 



probability oc exp 



(1) 



(A£)^ is an invariant, positive definite quadratic form which in the pair of 
independent thermodynamic parameters and may be written as 



2\2 



(2) 



where Ax" = {x°' — x^) {a = 1, 2) denotes the difference between the thermo- 
dynamic parameters x° of the subsystem and their values Xq corresponding 
to (A£)^ = 0. The thermodynamic metric elements ga/3 are evaluated in the 
state x" = Xq, and are tabulated below in Table 1 for several coordinate 
systems. 

The thermodynamic Riemannian curvature scalar (in the sign convention 
of Weinberg [8J) may be written a^ [3] 
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dx"^ \y/g dx^ y/g dx^ guy/g dx^ 



(3) 



where 



g = 911922 



^12- 



(4) 



^This equation is in footnote 53 of Ref. j3J. There is a small typographical error in Ref. 
[3] as the term dgi2/dx'^ should be dgi2/dx^. The equation is given correctly in {T,p) 
coordinates in Ref. [7]. 
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A{T, N,V) = U-TS 
n{T,ij, V) = U -TS - fiN 
G{T,p,N) = U-TS + pV 




dU^' dUdN' dm 
dSdN' dm\ 

d^VL d^Q d^VL 



d^G d^G d^G' 



keTV [dT^' dTdp' dp"^ 



Table 1: Thermodynamic potentials and the thermodynamic metric elements 
in five coordinate systems. Fluctuations for open subsystems take place at 
constant V, coinciding nicely with fixed V in the derivatives in the first 
four coordinate systems. In {T,p} coordinates, however, is held fixed 
on differentiating G{T,p, N). Nevertheless, the form of the {T,p} metric 
elements was picked to refiect fiuctuations at constant V [?]. Since R is 
an invariant, the choice of coordinates for calculating it is purely one of 
convenience. 
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R is an intensive thermodynamic variable with units of volume per molecule. 
Although the thermodynamic metric elements change their form on trans- 
forming coordinates, the value of R for a given thermodynamic state does 
not change since it is an invariant, by the rules of Riemannian geometry. For 
calculating R, the choice of coordinates is one purely of convenience. 

To conclude this section, I point out that Weinhold ^ originated thermo- 
dynamic energy metrics in the form of inner products based on the Hessian 
of the internal energy. The positive-definite nature of these inner products 
represents the second law of thermodynamics. But Weinhold's geometry 
lacks a true Riemannian metric structure since it has no underlying physical 
notion of distance, such as is offered by the fluctuation motivated entropy 
metric ^U\. An entropy metric was also used by Andresen, Salamon, and 
Berry [11] as a measure of the dissipated availability in finite-time thermo- 
dynamics. There have also been numerous calculations of R for black hole 
thermodynamics; see Aman et al. [12] for review. 



3 PHYSICAL INTERPRETATION OF R 

In this section, I present the basic thermodynamic curvature themes occur- 
ring in fluids. The discussion here expands the physical interpretation of R 
over what has been attempted previously. 



3.1 \R\ oc in the asymptotic critical region 

For the single component ideal gas, R = whether the ideal gas is monatomic 
or molecular. This finding originally motivated the hypothesis that R mea- 
sures intermolecular interactions [10] . The calculation of R in model systems 
quickly revealed a proportionality between \R\ and the correlation volume: 
\R\ oc with ^ the correlation length and d the spatial dimensionality, par- 
ticularly near the critical point where ^ diverges [31 UHl HSl HI]- oc in 
the asymptotic critical region is consistent with all the fluid data examined 
in this paper. 

To evaluate the dimensionless proportionality constant between \R\ and 
model calculations in which both R and ^ are evaluated are required. 
Several calculations of this type have been carried out in critical regions with 
large ^: 1) four pure fluids near the critical point [TO], 2) the one-dimensional 
ferromagnetic Ising model [I5] , and 3) the one-dimensional Takahashi gas [7] . 
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In all these cases, the same proportionality constant was obtained: 



R 



(5) 
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i? < in each case. Based on these limited calculations, the proportionality 
constant 1/2 would appear to be universal. I am not aware of any proof 
of this universality, and this is a matter which would seem to deserve some 
future attention. 

Following this start based on model calculations, the connection \R\ oc ^'^ 
was later confirmed on general grounds with a covariant thermodynamic fluc- 
tuation theory first developed by Ruppeiner [161 E] and completed by Diosi 
and Lukacs [1_8J who explicitly added the conservation laws. The idea is that 
at a large volume V the Gaussian thermodynamic fluctuation theory Eq. ([T]) 
works very well, but with decreasing V the fluctuating subsystem eventu- 
ally samples a correlated environment, which Eq. ([T| cannot model. The 
covariant thermodynamic fluctuation theory predicts that Gaussian fluctua- 
tion theory ceases to work at a volume [31 H] 



In the asymptotic critical region, V is physically interpreted as being roughly 
Clearly, this interpretation is consistent with the result Eq. (|5| based 
on direct calculation. For calculating the volumes of organized mesoscopic 
fluctuating structures below, I will use Eq. ^ since its derivation is based 
on general arguments. 

A pictorial depiction of the meaning of the correlation length, due to 
Widom [19], is given in Figure 2. Spontaneous density fluctuations cause 
the local density p (f^ at a point r in a single phase fluid to deviate from the 
overall density po in some complex, time dependent manner. Mathematically, 
p (r) = Po corresponds to an intricate contour surface separating two sides 
with local mean densities p (f) > po and p (r) < pg. A straight line through 
the fluid intersects this surface at points spaced an average distance ^ apart. 
^ is generally small in a disorganized system like an ideal gas, but diverges at 
the critical point for real fluids. This figure also shows a "droplet" of linear 
dimension ^ which offers schematic depictions of large spatially organized 
density fluctuations. Such a schematic droplet near the critical point is shown 
in Figure 5a. 



V 



R 



(6) 
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surface 



Figure 2: A surface in three dimensions on which the local density p {r) 
equals the overall density po- Also shown is an arbitrary line intersecting the 
surface at the dotted points, and a "droplet" of linear dimension ^ equal to 
the mean distance between those intersections. 

3.2 The sign of R 

A tabulation of results in a number of systems ^ suggests that the sign of R 
reveals the basic character of the intermolecular interactions; R is negative 
for thermodynamic states where attractive interactions dominate, and R is 
positive for states where repulsive interactions dominate. The clearest case 
is offered by the canonical examples of the Bose and Fermi ideal gasses, 
where quantum statistics causes atoms to either bunch closer together or 
further apart than in the corresponding classical ideal gas, mimicking the 
effect of attractive and repulsive interactions. Janyszek and Mrugala |2Cy, and 
Oshima, Obata, and Hara [21], showed that (in Weinberg's sign convention 
[H]) -R is always negative for the Bose gas and always positive for the Fermi 
gas. Mirza and Mohammadzadeh [22] worked out the ideal g-deformed boson 
and fermion gasses, which also show the appropriate sign. Brody and Ritz 
[13] examined the sign of R in finite Ising models as a function of the system 
size. 

Fluids are more complicated than ideal quantum gasses since the fluid 
intermolecular interaction potential typically has two parts, a repulsive part 
at short range and an attractive part at long range; see Figure 3. Sorting 
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Figure 3: Fluid intermolecular interaction potential 0(r) showing a repulsive 
force at short range, an attractive force at long range, and a potential well 
where molecules in the compact liquid or solid phases typically reside. 

out which part dominates in a given thermodynamic state can be difficult. 
One statistical mechanical guide is offered by the pair correlation function 



where I have assumed that the position vector r, with magnitude r, starts 
on a molecule within the bulk fluid. Fisher and Widom [23] argued that 
attractive interactions dominate if the long-range decay of G{r) is monotonic, 
and repulsive interactions dominate if this decay is oscillatory. These authors 
P3| were the flrst to calculate the Fisher- Widom curve along which the long- 
range decay of G{r) changes from monotonic to oscillatory. 

It is natural to ask whether the Fisher- Widom curve might coincide with 
the curve R = 0, since both these curves mark a transition from attractive 
to repulsive interactions. Near the liquid- vapor critical point, the long-range 
decay of G{r) is monotonic. The basic lore is that attractive interactions 
dominate, and organize fluctuations over long distances. One would thus 
expect negative R in the asymptotic critical region, and this is found in 
all the fluids examined here. As we move along the coexistence curve from 
the critical point towards the triple point, cases with R switching sign from 
negative to positive are found with the NIST Chemistry WebBook [5j in both 




(7) 
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the liquid and the vapor phases. However, convincing results concerning the 
long-range decay of G{r) are harder to come by. The possible correspondence 
between the R = curve and the Fisher- Widom curve is under investigation. 

3.3 Commensurate Ws along the coexistence curve 

Key in this research is the argument that the i?'s in the coexisting liquid 
and vapor phases are equal to each other in the asymptotic critical region. 
The structural underpinnings for this idea originated with Sahay, Sarkar, and 
Sengupta who calculated R for the van der Waals model, and observed 
that for isotherms with temperatures less than (but not too far from) the 
critical temperature, there must be an i?-crossing point where the curvatures 
in the liquid and vapor phases equal each other. This observation, augmented 
with a physical argument by Widom \19l, led to the proposal by Ruppeiner 
et al. |21] that this i?-crossing point coincides with the coexistence curve. In 
the Appendix, I offer a proof of this commensurate R theorem in the context 
of the currently accepted asymptotic scaling description of fluid criticality. 

The physical argument [23] for commensurate i?'s envisions approaching 
some point on the coexistence curve from either the liquid phase or the vapor 
phase; see Figure 4, where v is the molar volume. Widom [19] proposed that 
the correlation length in either of these bulk phases equals the thickness of 
the interface between the incipient coexisting phases. Since this interface 
thickness is the same at some point on the coexistence curve, no matter 
from which direction we approach it, the correlation lengths, and hence the 
i?'s, should be equal in the coexisting liquid and vapor phases. This is the 
case even though the critical point joining the liquid and vapor phases is a 
singular point, and even though the coexisting liquid and vapor densities can 
be quite different from each other. I add that Evans et al. ^6J showed with 
density functional theory using a short range intermolecular fluid potential 
that the character of the density decay, including the correlation length, at 
the interface matches that in the bulk, lending even further support to the 
Widom argument. 

This idea offers a practical new way of dealing with first-order phase tran- 
sitions. As is well known, ^ can be difficult to calculate or to measure but R 
follows readily from thermodynamic properties. In the van der Waals model, 
the commensurate R theorem was used in place of the physically problematic 
Maxwell equal area construction to derive the liquid-vapor coexistence curve 
[24j . May and Mausbach [27] followed up by calculating the liquid- vapor co- 
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Figure 4: The liquid-vapor coexistence curve and the critical point. Accord- 
ing to the commensurate R theorem, at a given temperature T, the i?'s in the 
two coexisting phases are equal, even though their molar volumes v might 
differ considerably. 

existence curve for Lennard- Jones computer simulation data. In both cases, 
strong results were obtained in this otherwise difficult problem. 

Widom [19] also argued that if we have a majority phase of vapor on the 
verge of a first-order phase transition, then the fluctuating droplets within 
it have the density of the liquid phase which is to be made. Likewise for 
the majority liquid phase preparing to make the vapor phase. Figure 5c 
illustrates both this idea and the commensurate R theorem. 

3.4 Low \R\ limit 

Outside the asymptotic critical region, a limiting theme enters the picture 
|24j . With decreasing \R\, the correlation volume could become less than 
the molecular volume v, particularly in the vapor phase where v typically 
becomes very large as the triple point is approached. In this low \R\ limit, 
with \R\/v -C 1, we run short of statistics to do thermodynamics with at 
volume scales of \R\, and attempts at physical interpretation of R might 
appear unreasonable. Although such physical interpretations are attempted 
here anyway, the possibihty that our use of R might be stretched beyond its 
limits must always be considered. 
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Figure 5: Schematic figures illustrating several mesoscopic fluid concepts: 
a) a loose organization of molecules with volume \R\ pulled together by the 
attractive part of the intermolecular interactions (i? < 0), b) a compact 
cluster of molecules of volume \R\, pulled together by the attractive part of 
the intermolecular interactions, but prevented from collapse by the repulsive 
part of the intermolecular interactions {R> 0), c) a fluid in two phases near 
the critical point, with the bottom half the liquid phase containing vapor 
droplets with volume \Ri\, under a coexisting commensurate vapor phase 
containing liquid droplets with the same volume Ry = Ri as those in the 
liquid, d) hquid and vapor phases with incommensurate droplet sizes, e) a 
disorganized compact liquid phase held together by attractive intermolecular 
interactions {R < 0), and f) an organized compact solid phase held up by 
the repulsive part of the intermolecular interactions {R> 0). 

3.5 Incommensurate Ws along the coexistence curve 

The commensurate R picture associates liquid-vapor phase coexistence with 
mesoscopic structures of naturally the same spatial size ^ in the two phases. 
I conjecture that if these structures were of dissimilar sizes, as shown in 
Fig. 5d, and as is typically the case outside the asymptotic critical region. 
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then it might be difficult for one phase to form the other, corresponding to 
metastabihty. 

Particularly difficult to treat theoretically is the onset of boiling as a 
liquid is warmed. The classical homogeneous nucleation theory |28] is shown 
in Figure 6. A bulk liquid at temperature T contains a vapor bubble of radius 
r and pressure attempting to expand against the sum of the liquid pressure 
pi and the pressure caused by the bubble's surface tension S {S > Qi). It is 
generally assumed in this theory that all thermodynamic properties relate in 
the same way as in the thermodynamic limit, questionable at length scales 
approaching the order of intermolecular distances. 

For the bubble to grow, and form a macroscopic vapor phase, we require 

[28] 

25' 

Pv>Pi^ , (8) 

r 

where pi is set by some external constraint, such as atmospheric pressure. 
As the bulk liquid is warmed to the coexistence curve, PviT) — ^ Pi- This 
limit cannot satisfy Eq. ([s]), and, for the liquid to boil, we must superheat 
by warming to a T higher than the phase transition point, with a larger p^. 

The homogeneous nucleation theory clearly cannot work for tiny bubbles, 
since 2S/r increases without limit as r — )■ 0. Therefore, we usually assume 
that the bubbles are initially larger than some critical radius r = hy 
forming on nucleation sites in the bulk liquid or on the container walls. S 
generally decreases to zero at the critical point [19j, and so we would not 
expect much superheating in the asymptotic critical region. 

Temperatures below which superheating might become significant have 
been associated with corresponding isotherms having negative pressures in 
the metastable liquid regime. In such the liquid is able to withstand a 

pulling stress in the related cavitation problem [28]. In simple fluid models, 
such negative pressures typically take place at T < 0.9Tc [28] (for van der 
Waals, T < 27X^/32 [2^). 

The homogeneous nucleation theory of boiling is conceptually challenged 
since, as becomes clear in section 4, for T ~ 0.9 Tc, typically |-R| ~ Inm^ in 
the liquid phase. Such a volume contains roughly one vapor molecule, and 
it is hard to see how the macroscopic rules of thermodynamics assumed in 
this bubble model could possibly hold. Zeng and Oxtoby [30] gave a critical 
analysis and concluded that " although the standard classical theory has been 
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Figure 6: Homogeneous nucleation theory picture of boiling. A bubble of 
vapor with radius r and pressure attempts to expand inside a liquid at 
temperature T. To expand, the bubble must overcome the sum of the bulk 
liquid pressure pi and the pressure 2S/r contributed by the surface tension 
S. 

quite successful in predicting droplet formation for many materials, it should 
fail completely in the process of bubble formation." 

The suggestion in this paper that incommensurate i?'s form a barrier to 
the onset of the liquid-vapor phase transition puts the focus on fluctuating 
mesoscopic structures, and without assuming anything about their particular 
character. 

3.6 R in the compact liquid phase 

The NIST Chemistry WebBook data presented in section 4 indicates that 
the compact liquid phase near the triple point has small \R\, on the order of 
the volume of a molecule. The sign of R is mostly negative; the few instances 
with positive R coincide with having a curve with R = intersect the liquid 
phase. 

The only continuum fluid model where results for R in the compact liquid 
regime have been reported is the one-dimensional Takahashi gas [7J. These 
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results support small of the character found here. I add to this calculation 
with a simple hard-sphere liquid model having Helmholtz free energy 

A(T,iV,y) = iVA;BT In (^-^^ +iVe(T)- iV^BT In (^1-^^ , (9) 

where p is a constant with units of density, e(T) is a function of T with 
negative second derivative, and the constant b gives the hard-packing limit. 
The particle density is p = N/V, the pressure p = pkBT/{l — bp), and the 
heat capacity at constant volume Cy = —NTe"{T). By Table 1 and Eq. (|3|, 

R = -b{l-bp). (10) 

Clearly, — 6 < i? < in the regime of physical densities < p < 1/6. 

In the covariant version of the theory, \R\ sets the lower limit of applicabil- 
ity of thermodynamic fluctuation theory [16]. Since for the hard-sphere gas, 
thermodynamic fluctuation theory is expected to work down to length scales 
on the order of the volume of a molecule 6, the result \R\ ~ 6 is expected. 

This somewhat primitive model gives little insight into the sign of i?, 
which in this model is uniformly negative. To discuss the sign of i?, I return 
to real fluids. Near the critical point R is always negative in the liquid state. 
As 1 show with the real fluid data in section 4, as we cool the liquid state 
towards the triple point, \R\ typically decreases until it takes on a value of the 
order of the volume of a molecule. Figure 5e pictures such a compact liquid 
state with i? < 0, loosely held together by attractive interactions. We now 
have the question: as we cool further, could this state organize itself into a 
"solidlike" structure where hard-sphere repulsive interactions dominate and 
where R changes sign to i? > 0, as shown schematically in Figure 5f? As 
seen in section 4, such a sign change is indeed present in a few fluids. 1 add 
that Widom [31] contrasted the liquid behavior at the critical point and the 
triple point, and made the point that the former is dominated by long-range 
attractive interactions, and the later by short-range repulsive interactions. 

3.7 in the vapor phase 

For the vapor phase R is negative in the asymptotic critical region. As we 
cool along the coexistence curve from the critical point to the triple point, we 
expect to decrease at first, and eventually end up negative near the triple 
point. Negative R is to be expected in any vapor regime where the molar 
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volume is large, the molecules far apart, and the long-range attractive part of 
the intermolecular potential dominant over the repulsive part. Negative i?'s 
near the triple point are indeed seen in all the vapors considered in section 
4. 

The magnitude of R in the vapor phase at the triple point is, however, 
harder to interpret physically than its sign. The fluids in section 4 show 
that either the vapor R settles down to about — Inm^ or decreases to large 
negative values. The later case generally occurs in vapors with very large 
triple point molar volumes. In either of these scenarios, the triple point 
\R\/v ^ 1, and we are clearly below the length scale set by the low \R\ 
limit. It would thus be easy to dismiss these vapor phase triple point results 
as being devoid of physical significance. However, these results do not stand 
alone since qualitatively similar results were obtained for the one-dimensional 
Takahashi gas at large molecular volumes [7]. But I attempt no physical 
interpretation here. 

As we cool from the critical point, an interesting new feature presents 
itself in the vapor phases of a few fluids, such as water. In water near the 
critical point, but outside the asymptotic critical region, R crosses zero twice, 
with a peak at a positive value i? ~ +1 nm^. I associate the resulting regime 
of i? > with the cluster forming scenario pictured in Figure 5b. 



4 NIST CHEMISTRY WEB BOOK RESULTS 

In this section I report results of calculations of R using fluid data from the 
NIST Chemistry WebBook [3] . The fluids I examined are made of molecules 
of four types: 1) monatomic molecules: Hehum [32], Neon [33], Argon 
Krypton [35], Xenon [35], and Methane [36], with Methane placed in this cat- 
egory because it is quasi-spherical [37], 2) linear diatomic molecules: Normal 
Hydrogen [38], Nitrogen [39], and Oxygen [ID], 3) other linear molecules: 
Carbon Dioxide [H] and Carbon Monoxide [35], and 4) more complicated 
molecules: Water [H], Methanol [43j, and Hydrogen Sulfide [33] . 

My calculations represent a broad based approach to significant themes 
for thermodynamic curvature in fluids. Such a broad approach certainly fits 
the spirit of the NIST Chemistry WebBook, which is built on correlations, 
averaging, and extrapolations using numerous data sets in many fluids. The 
NIST Chemistry WebBook is, however, not optimal in all circumstances. For 
example, near the critical point it is not based on the scaled critical equations 
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of state, and I will thus attempt no detailed critical point analysis, such as 
the evaluation of the critical point exponent of R. But, fortunately, such 
analysis is not necessary in a first approach, given the critical point theorems 
presented in the Appendix. 

The basis of the calculation of R is the Helmholtz free energy per volume 



/(T,p) = 



A{T,N,V) 
V 



(11) 



whose numerical values may be looked up in the NIST Chemistry WebBook 
[3] 1^ The NIST Chemistry WebBook is based on fitting experimental data to 
smooth multiparameter formulas, yielding precise numbers very well suited to 
calculating numerical derivatives with finite difference formulas [H]. In (T, p) 
coordinates. Table 1 yields immediately the thermodynamic line element 
Ae = QttAT^ + QpAp"^ with 



{9TT,9pp} 



1) 



(12) 



In these coordinates, Eq. ([s]) becomes 



R 



d 



^ [dT \^ dT 



^ dp dp ) 



with 



g = 9tt Qpp- 



(13) 



(14) 



Hydrogen offers a nice first case, with its thermodynamic properties re- 
cently reanalyzed [5B] • Figure 7a shows R for Hydrogen along the coexistence 
curve in both the liquid and the vapor phases, from the critical point to the 
triple point. A log-log plot of R versus (Tc — T)/Tc shows a power law diver- 
gence {R — i- — oo) at the critical point with slope near 2 in both the liquid 
and vapor phases [21]. Such a divergence is seen at least approximately for 



all the fluids examined in this paper, and is expected from Eq. (28). R for 



the liquid and vapor phases agree with each other to better than 1% in the 
temperature range 0.96 < T/Tc < 1, consistent with the commensurate R 
theorem. By contrast, at T/Tc = 0.96, the molar densities of the coexisting 
liquid and vapor phases differ from each other by a factor of ~ 3. 



^ Since the data tabulations are in quantities per mole, a multiplication by p is required 
to convert to quantities per volume. 



17 



|a) Norinal Hydfoganl 



7"c = 33.19K 



Tf - 13.96K 18.8K 



20 25 
T(K) 




-10 



b) Helium I 



F;, = 2.177K 



2.0 2.5 3.0 



Tc = 5,195K 



3,80K 



LiqoltfR 

Vapor R 

• R = 



3.5 4.0 4.5 5.0 



1-1 



-s ■ 



I c) Argon | 



T( = 83.81 K 




Tr = 190. 56K 




80 90 100 110 120 130 140 150 80 100 120 140 160 180 200 



r(K) 



Tj= 54.361K 



a) Oxygan | 




-10 



| f) Water 



T, = 273. 16K 




Liquid R 

Vapor R 

• K-0 



60 80 100 120 14D 160 250 300 350 40O 450 500 550 600 650 

r(K) T(K) 



Figure 7: R for six representative fluids along their coexistence curves: a) 
Hydrogen, b) Helium, c) Argon, d) Methane, e) Oxygen, and f) Water. The 
data go from the critical point, with temperature Tc, to the triple point, 
with temperature Tt (except Helium which starts at the lambda point with 
temperature Tx). Each case shows R diverging to —oo at the critical point, 
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with the it!'s more or less commensurate in the two phases in the asymptotic 
critical region. 
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Fluid 


' t 


Vl 




Ri 


hi 


Rilvi 


7?^, / V^, 
' / 






01 fil 1 

U.U-LU ± X 




-0 lf^f^4 

-U. i.yjO'-t 


-0 XTQ 




-0 0^0 




00.0 ± 






-0 09^0 


-0 QR*^ 


-0 


-0 OfiO 


Krypton 


115.78 


0.03425 


12.742 


-0.0286 


-1.697 


-0.502 


-0.080 


Xenon 


161.41 


0.04426 


15.962 


-0.0373 


-2.236 


-0.508 


-0.084 


Methane 


90.69 


0.03553 


63.925 


-0.0113 


-1.049 


-0.192 


-0.010 


Hydrogen 


13.96 


0.02618 


15.501 


+0.0083 


-0.5407 


+0.190 


-0.021 


Nitrogen 


63.15 


0.03230 


41.483 


-0.0162 


-0.7671 


-0.303 


-0.011 


Oxygen 


54.36 


0.02450 


3,081.0 


+0.0012 


-141.99 


+0.029 


-0.028 


CO2 


216.59 


0.03735 


3.1972 


-0.0694 


-0.6276 


-1.119 


-0.118 


CO 


68.16 


0.03297 


36.091 


-0.0177 


-0.9041 


-0.322 


-0.015 


Water 


273.16 


0.01802 


3,710.9 


-0.0018 


-114.01 


-0.059 


-0.019 


Methanol 


175.61 


0.03542 


7,834,000 


-0.0130 


-689,430 


-0.221 


-0.053 


H2S 


187.70 


0.03434 


66.557 


-0.0155 


-1.2778 


-0.272 


-0.012 



Tabic 2: Triple point table showing the triple point temperature Tt in Kelvin, 
the molar volumes in the liquid and vapor phases, vi and v^, in liters per mole, 
the liquid Ri and the vapor R^, in nm'^ per molecule, and Ri/vi and R^/vy 
dimensionless. 

For Hydrogen \R\ in the liquid phase decreases on cooling towards the 
triple point, with R changing sign and remaining positive below T = 18. 8K. 
At the triple point, R = +8.3 x 10^'^nm"^; see Table 2. A corresponding sphere 
with volume |i?|/6 would have radius 0.69A, of the same general order as the 
van der Waals radius of a Hydrogen atom 1.2A. An organized structure with 
this size is in accord with the expectations from section 3.6 for the compact 
liquid phase. As is clear from Table 2, such liquid phase triple point values 
of \R\ are typical for all the fluids. Most of the liquid phase triple point i?'s 
are negative, characteristic of compact disorganized structures, as in Fig. 5e, 
but there are three fluids with positive R, characteristic of the order in the 
solidlike states, as in Fig. 5f. 

In the vapor phase for Hydrogen, R is uniformly negative. This is rea- 
sonable, since in the asymptotic critical regime R is negative, and as we cool 
towards the triple point the molecules become ever more widely spaced, and 
the attractive part of the interactions will thus increasingly gain in signifi- 
cance over the repulsive part. Negative R in the vapor phase is characteristic 
of all the fluids I looked at except Water and Methanol. (Neon vapor has 
positive R near the triple point, but this may be spurious.) 

R for Helium, shown in Fig. 7b, looks qualitatively similar to that for 
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Hydrogen, though the coexistence curve for Hehum ends not at a triple point 
but at a lambda point at temperature Tx = 2.177K. Again, in the asymptotic 
critical region, we see strong agreement with the commensurate R rule. In 
the vapor phase, R is uniformly negative, but in the liquid phase R changes 
sign and become positive below T = 3.80K. 

R for Argon, shown in Fig. 7c, has the commensurate R theorem not as 
conspicuously satisfied as it was for Hydrogen and Helium. Whether this re- 
flects the real behavior, or shortcomings in the fitting formula in this regime 
is not clear. For Argon, both the liquid and the vapor i?'s are uniformly 
negative along the entire coexistence curve. The basic pattern for Argon 
is repeated for Krypton, Xenon, Methane (see Fig. 7d), Nitrogen, Carbon 
Dioxide, Carbon Monoxide, and Hydrogen Sulfide. Neon shows similar be- 
havior as well, except for a regime of small positive R in the liquid phase 
near the triple point. Neon terminates in an anomalously large liquid \R\/v, 
as shown in Table 2. 

R for Oxygen, shown in Fig. 7e, has the liquid phase terminate with a 
positive R at the triple point. R changes sign at T = 57.9K, with \R\/v 
at the triple point anomalously small; see Table 2. Whether or not this is 
a real effect, or reflects a problem with the data fit, is unclear. The vapor 
phase shows a new feature, with R decreasing abruptly on approaching the 
triple point. As is indicated in Table 2, such large vapor phase \R\ values are 
typically associated with corresponding large increases in the vapor molar 
volume. Water shows a similar feature, which may be related to metastabil- 
ity. 

R for Water, shown in Fig. 7f, shows typical liquid phase behavior. How- 
ever, the vapor phase shows a new feature which may be physically signifi- 
cant, a positive peak for R at T = 616K, with R = 0.341 nm^. (For this state 
R/v = 1.1). I associate this feature with cluster formation of the type shown 
in Fig. 5b, where water molecules are in a condensed solidlike state held up 
by the repulsive part of the intermolecular interactions. It is straightforward 
to estimate the average number of water molecules in a cluster at the peak 
R. By Eq. ([6]), the cluster volume is roughly \R\/Q = 0.0568 nm^. Assume 
that each water molecule occupies a spherical volume with van der Waals ra- 
dius: (47r/3)(0.17nm)^ = 0.021nm^, leading to a number of molecules in the 
cluster of 0.0568/0.021 ~ 3 molecules. Computer simulations by Johansson 
et al. [l5] have found a significant number of water clusters in the coexisting 
vapor corresponding roughly to this peak R. These clusters were found to be 
dominated by dimers, though clusters as large as 7 molecules were identified. 
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Methanol shows a similar feature of positive R in the vapor phase near 
the critical point, but the quality of the NIST fit may not be as good as for 
Water. If this connection between positive R and physical clusters is correct, 
then calculating R offers a method of identifying fluids and thermodynamic 
regimes where we might find clusters. This idea will be pursued elsewhere. 

For Water, there is a point at T = 529K where the vapor curve for R 
crosses the liquid curve, and a rapid divergence between the values of R 
in the two phases occurs as we cool to the triple point. If we interpret 
these large R differences in terms of metastability, then we clearly see that, 
as we warm from the triple point, metastability should end at T = 529K. 
Brennen [2S| gives the maximum temperature of superheating experimentally 
observed in Water as about 550K, with corresponding van der Waals value 
T = 27/32Tc = 546K. 

5 CONCLUSION 

In conclusion, this overview of the thermodynamic curvature R on the liquid- 
vapor coexistence curve of fourteen pure fluids has confirmed some old ex- 
pectations, but also brought forth some new features. The broad trend is 
negative values of R, as the attractive part of the intermolecular interactions 
usually dominates for pure fluids, with the exception of compact situations 
where the molecules are in close contact. In the asymptotic critical region, 
|-R| oc with the liquid and vapor values of R in the coexisting phases equal 
by theorem. On cooling towards the triple point, \R\ typically decreases until 
it is of the order of molecular sizes in the liquid phase, but less well defined 
in the vapor phase. 

Interesting new features are of three types. The first is incommensurate 
i?'s outside the asymptotic critical regime, with fluids of largely different i?'s 
in the liquid and vapor phases. I associate this with metastability, where 
the majority phase has a difficult time making the minority phase because 
the mesoscopic fluctuating structures in the phases are of different sizes. 
Although I present no detailed experimental examination supporting this 
conjecture, it would appear to be at least logical. Second, I associate isolated 
regimes of positive R in the vapor phase near the critical point with the 
formation of solid clusters. Third, fluids with R changing sign through R = 
in the compact liquid phase may be associated with the Fisher- Widom curve 
marking the transition in the long-range decay of the correlation function 
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from a monotonic fluidlike character to an oscillatory solidlike character. 

Generally, further study of R might provide insight not only into the 
nature of thermodynamic curvature, but add another set of criteria to use in 
judging the quality of fits of fluid data to fitting models. 



I thank Peter Mausbach, Helge-Otmar May, Horst Meyer, Anurag Sahay, 
Tapobrata Sarkar, Guatam Sengupta, and Steven Shipman for useful com- 
munications. Very special thanks to Eric Lemmon for programming R into 
REFPROP 9.01 (/3 test version). 



7 APPENDIX: COMMENSURATE R THE- 
OREM 



Significant aspects of this paper hinge on the behavior of R in the asymptotic 
critical region. Below, I work in the context of the scaled form of the free 
energy, with the Rehr and Mermin |16] mixing of scaling variables to deal 
with asymmetry. This constitutes the " currently accepted asymptotic scaling 
description of fluid criticality" 07]. In this Appendix, I prove that R has the 
same critical exponent as along the coexistence curve and its analytic 
extension into the supercritical regime. I also prove the commensurate R 
theorem that R in the coexisting liquid and vapor phases are equal. 

My proofs assume that along the coexistence curve n = /i(T) is analytic 
approaching the critical point from below. This assumption has been ques- 
tioned in the context of the Yang- Yang anomaly where the second derivative 
/i"(T) diverges at the critical point. There is some experimental evidence 
[I8] for this point of view which requires "complete scaling" to address the- 
oretically. In this scenario, my proofs would require revisions. 

Write the pressure as [i6] 



where po{T, /i) is the regular analytic part, and the term containing the func- 
tion Y±{z) is the singular part. 
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p(T,/x)=po(T,/i) + |rrF±(z) 



(15) 




(16) 
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(a, b) are two critical exponents, 



r = (T-Te)/T, + ci(/i-/ie)ABT„ (17) 

C={fl- lXc)/kBTc + C2(T - Te)/T„ (18) 

{Tc, fic) are the values of (T, /i) at the critical point, and ci and C2 are two 
constants. I assume |ciC2| < 1, the case if the fluid is not too antisymmetric. 
The function Y±{z) has no explicit dependence on ci and C2, and has two 
branches (±) depending on the sign of r. These branches join smoothly 
along the curve r = 0, except at the critical point {t,(} = {0,0}. The 
critical exponents a and b are related to the standard critical exponents a 
and /3 by a = 2 - a and 6 = UHl HQ]. 

It is assumed that Y±{z) is an even function of z. This assumption, and 
the goal of modeling a first-order phase transition terminating in a second- 
order critical point, requires [46]: 

n(o) = 0, (19) 

yI^\0) = 0, (20) 

r4o+) = vi(o_)7^o, (21) 
r(o+) = -r(o_)^o, (22) 



r'(o+) = r'(o_) ^ 0, (23) 



and 



y1^\0+) = -Y^^\0_) 0, (24) 



where the subscripts -|- and — on refer to the sign of z as either ^ — 0+ or 

C^o_. 

For T < 0, C = 0± represents the two phases on the coexistence curve 



because, for given (T, /x), we have p(r, 0+) = p(r, 0_), by Eqs. (15) and 



(21). This continuity of (T, p, //), and the corresponding discontinuities of 



the density p = {dp/dp)T and entropy per volume s = (dp/dT)^ resulting 



from Eqs. (15) and (22), are the conditions for a first-order phase transition. 
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Note that setting ^ = O-t in Eq. (18) yields (/i — Hc)/kBTc = — and 



r = (1 — ciC2)t in Eq. (17), where the reduced temperature 

T - T 

t= " 



(25) 



Thus, on the coexistence curve, t is negative and {n — /ksTc has the same 
sign as C2. 

In (T, /i) coordinates the thermodynamic hne element can be read off 
from Table 1: 



Qf2 



AT' + 2 J^ATA/i + I^A/i^ 



(26) 



where I have set Q = —pV, and T = Tc in the prefactor since we are in the 
asymptotic critical region. It is straightforward to show that the term con- 



taining Y±{z) in Eq. ( 15 ) causes each second derivative in the thermodynamic 



line element Eq. (26) to diverge on approaching the critical point, assuming 
critical exponent values in the vicinity of the pure fluid ones (a ~ 1.9, b ~ 1.6 
poj). Therefore, since the regular part of the pressure po{T, fi) produces only 
finite second derivatives in the metric elements, we can drop it in calculating 
R in the asymptotic critical region. 

For r > 0, and along the line ( = representing the analytic continuation 
of the coexistence curve, Eq. ([s]) yields. 



kBT,ib-l){2b-a)il-c,C2] 
a(a -l)r+(0) 



(27) 



demonstrating that R+{t, 0) diverges with exponent a = 2 — a. This exponent 
is the same as that of [I9j. To assure a positive heat capacity, we must 
have l+(0) > 0, yielding R+{t, 0) < for fluids. This is found in every case 
considered in this paper. 

Along the coexistence curve we have exactly 



-A;bT,(6- 1)(1 
x[-(a-6)2(H 



6-a)r (0±)2r'(0d 



i?_(t,0±) 



+ 2a(a- l)(26-a)yi(0±)r'(0±)2 

+ a(a- l)(a-6)F_(0±)F:(0±)F?^(0±)] 

2 [{a - 6)2r:(0±)2 - a{a - l)Y_{0±)Y!!{0±)f 



(28) 
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Every product of Y_ {0±) and its derivatives has the same value in either phase 



by Eqs. (21)-(24), estabhshing immediately the commensurate R theorem: 

i?_(t,0_) = i?_(t,0+). (29) 
Note also that i?„ diverges with exponent a = 2 — a, the same as that 



of The presence of the term yP{0±) in Eq. (|28|), with sign unset by 



thermodynamic stability, means that Eq. (28) does not clearly set the sign 
of R- {t, 0±). But it was found to be negative in the asymptotic critical region 
in all fluids examined in this paper. 

Other than T, p, /i, and R, most thermodynamic functions will not be 
equal in the coexisting phases. For example, the metric elements ga/s, related 
to heat capacities and compressibilities, consist of a mixed sum of terms, 
some changing sign on switching phases and some not changing sign. 
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